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Abstract. We present two formulations of Fokker-Planck-Roscnbluth-type 
(FPR) equations for many-particle self-gravitating systems, with first order 
relativistic corrections in the post-Newtonian approach (1PN). The first starts 
from a covariant Fokker-Planck equation for a simple gas, introduced recently by 
G. Chacon- Acosta and G. Kremer [S]. The second derivation is based on the 
establishment of an 1PN-BBGKY hierarchy, developed systematically from the 
1PN microscopic law of force and using the Klimontovich-Duprec (KD) method. 
We close the hierarchy by the introduction of a two-point correlation function that 
describes adequately the relaxation process. This picture reveals an aspect that is 
not considered in the first formulation: the contribution of ternary correlation 
patterns to the diffusion coefficients, as a consequence of the nature of 1PN 
interaction. Both formulations can be considered as a generalization of the equation 
derived by Rezania and Sobouti in 2000 pQ , to stellar systems where the relativistic 
effects of gravitation play a significant role. 



PACS numbers: 05.20.Dd, 05.10.Gg, 04.40.-b, 04.25.Nx 
1. Introduction 

The collisionless Boltzmann equation in 1PN approximation for a self gravitating 
system, imbedded in an otherwise flat space-time, was derived by Rezania and 
Sobouti in 2000 pQ, finding some relevant solutions. The aim of the present paper is to 
go beyond the collisionless case and incorporate situations where encounters between 
particles play a significant role, and derive tractable kinetic equations describing the 
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evolution of this class of self gravitating systems. As it was pointed out by Kandrup, 
the astrophysical objects where the relativistic effects could play a significant role 
in the relaxation process are the galactic nucleus and, perhaps, the relativistic star 
clusters with age shorter than its relaxation time [TT]. There is solid observational 
evidence that there exists nuclei containing massive black holes formed due to the 
dynamical instability of relativistic systems of stars. In particular, the relativistic 
effects might be important in the case of a galactic nucleus decoupled from the 
remainder of the galaxy, which becomes even more dense (and relativistic) at a time 
scale of the order of its relaxation time [21 [3] . Since these systems are composed by 
stars with velocity much less than the speed of light, the post-Newtonian approach 
supplies the adequate tool to investigate their behavior. 

It is well known that the evolution of typical globular clusters, where it is 
assumed that relativistic corrections are unimportant, is described satisfactorily by 
the Fokker-Planck kinetic equation. In the local approximation, it is possible obtain 
explicit relations for the diffusion coefficients, in terms of Rosenbluth potentials, and 
such equation takes a very tractable form [I] (in this paper we call it as the FPR 
equation) : 

where we have used the sum convention = 1,2,3). Here A l N and B 1 ^ are the 
diffusion coefficients or Rosenbluth's coefficients (the subscript N remarks the fact 
that we are dealing with the classical Newtonian theory). Such equation describes 
the evolution of the distribution function /(x, v,t) of a test star, moving through 
a homogenous sea of field stars. They are described by a static space-independent 
distribution function, ^(v), which determines D i and D l K If each field star have the 
same mass m than the test star, the Rosenbluth coefficients take the form [U [5] 

A% = SvrGWlnA^- / d 3 v'^^-, (2) 
oviJ |v-v'| 

B% = 47rG 2 m 2 lnA-^— f dVtf(v')|v - v '|, (3) 

ovidvj J 

where In A is the Coulomb logarithm and G the gravitation constant. As it is well 
known, derivation of ©-© requires three fundamental facts: (i) To assume that the 
binary weak encounters dominate the system's relaxation process (6j[7]; (ii) to assume 
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that all of these encounters are local (i.e. with impact parameter much less than the 
system size); (iii) a detailed knowledge of the two-body microscopic dynamics. 

Perhaps, a natural way to obtain an 1PN version of (pQ), is to start 
from a covariant Fokker-Planck equation and then perform its post-Newtonian 
approximation. At present, it can be find different versions of this kind of kinteic 
relations, such as the one derived by Kandrup P, HO, E] based on a stochastic 
standpoint (i.e. from a covariant master equation). A more fundamental approach 
was introduced recently by G. Chacon-Acosta and G.M. Kremer [8j, who derived a 
covariant Fokker-Planck equation for a simple gas in the presence of gravitational field 
(see eq. fll])), starting from the general-relativistic Boltzmann equation. There are 
several reasons impelling us to adopt this scheme. Since such derivation is supported 
on the Boltzmann equation, the corresponding Fokker-Planck relation is consistent 
with situations where the system evolves towards an equilibrium sate, characterized 
by the Maxwell- Jutner distribution function |12j . On the other hand, such equation 
is valid for systems whose relaxation process is dominated by grazing collisions, which 
is the usual picture to model the behavior of self-gravitating systems (in this context 
it is more convenient to use the term encounters instead of collisions). More over, 
the diffusion coefficients characterizing (Tj0) are presented in a form that facilitates 
the implementation of the post-Newtonian approach. We shall show this in section 
El dedicated to obtain the 1PN approximation of (j4j). The result is an equation 
that is not manifestly covariant (this is a characteristic of the 1PN scheme), whose 
convective (l.h.s) and collision (r.h.s) terms can be split on a Newtonian and post- 
Newtonian contributions. 

The limitations of the above scheme can be viewed through the consideration of 
a more fundamental standpoint, i.e. a BBGKY formalism [THJ [TTl [THl [T9] , which will 
be the subject in later sections. By examining the system's dynamics at a microscopic 
level, one can note that a fundamental aspect of the 1PN interaction has not been 
contemplated in the formulation of (00). As we shall show in section [3], the 1PN force 
exerted on each particle has a ternary nature (see (l35l) . (!38|) ). This feature leads to 
the apparition of 3-order correlation patterns in the first BBGKY equation and not 
all of those vanish after the implementation of binary encounters assumption (see 
( l49l) . ( |52|) ). which is the basic statement employed in the formulation of collision 
terms of the Landau [TBI [HJ [15] or Fokker-Planck type. 

Despite of the apparition of third order correlation patterns originates a collision 
term that does not look like in a standard fashion (eq. ( 1521) ). it can be put in a 
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FPR form by choosing adequately the two point correlation function g 2 and setting 
#3 = 0. Here, two facts will play an important role: (i) The classical FPR equation 
can be derived from the first BBGKY equation by choosing a correlation function 
corresponding to a weakly coupled gas (WCG) in the hydro dynamical regime (it is 
not a surprising fact, since in the usual derivation of the FPR equation it is assumed 
that the weak encounters play a dominant role in the relaxation process); (ii) in the 
1PN approach the momentum conservation law for a system of point particles is 
the same as in Newtonian theory, and it holds if and only if each particle obeys the 
Newtonian equation of motion [22]. This means that the scattering process and, in 
consequence, the explicit form of g 2 in 1PN approximation are the same as in the 
Newtonian case. These considerations permit us to incorporate the hydrodynamical 
WCG correlation function in the post-Newtonian contribution of the collision term 
for its subsequent simplification. 



2. Derivation from a Covariant Fokker-Planck Equation 

We start by considering that the distribution function f(x' I ,p' 1 ) of a self-gravitating 
gas of particles with identical rest mass m, satisfies the covariant Fokker-Planck 
equation [8]: 

^-W^ = -^(/^) + ^(/^) (4) 

where p M = mU^ is the four-momentum, T° are the Christoffel's symbols, x M = 
(ct, x) is the set of configuration coordinates and D^, D^ v are the diffusion coefficients 
given by 

D» = [ /.Atfi^tftV^— , (5) 
J P*o 

D>» = [ UAp^Fo-dQV^g—, (6) 
J P*o 

where /* = /(x'Spj), Ap% = p'f - p£, g = det(g^), F = y/ (p%p a ) 2 - m 4 c 4 is the 
invariant flux, a is the differential cross section and dQ is an element of solid angle 
which characterizes the scattering process. The quantities p M and p^ denote the four- 
momentum of the test and field stars before the encounter, respectively, while 
and p'f represent their four-momentum after the encounter. 
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In order to perform the 1PN approximation of (jl]), we have to take into account 
the following facts. We need an expansion of (BJ-(l6l) up to order (w/c) 4 , where v 
is a typical Newtonian speed in the system and c is the speed of light. In this 
approximation g^ v is given in terms of the Newtonian field <3> and the post-newtonian 
fields if), £i (see |Appendix A[ ) : 

2oo= -l-2$/c 2 -2($ 2 + ^)/c 4 , (7) 
9oi = &/c 3 , (8) 
9ij = (1 - 2$/c 2 )cV (9) 

(Latin indices run from 1 to 3). On the other hand, the four- velocity W 1 is related 
to the classical velocity v % = dx l jdt through the equation 

U» = U°V fl , V = (l,v i /c), (10) 

and it is restricted by the relation 

giiV U»U" = -c\ (11) 

Another important fact is that in the 1PN approximation the momentum 
conservation for a system of point particles is satisfied if and only if each particle 
obeys the Newtonian equation of motion [22]. In other words: they obey the 
classical momentum conservation law. Immediately we note two implications that 
will play a fundamental role in the calculation of the diffusion coefficients. They were 
obtained in reference [8] by choosing the center-of-mass system corresponding to two 
colliding particles. Since in 1PN approximation they obey the classical momentum 
conservation law, it implies that AU® = and, in consequence, D° = D 00 = 0. 
On the other hand, also in virtue of the classical momentum conservation law, we 
must to consider adfl as the usual Newtonian scattering cross section. The above 
considerations enable us to rewrite (SJ)- (J6J) as 

M = -mV A ' )+ tm>mV* t) (12) 

where Cu is the Liouville's operator defined as 

C ° = u *t*-W u "w (13) 

and 

/d 3 U 

f*AuW(u?u a y - c^adn^—^, (u) 
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f*AU*AUW(U?U a ) 2 - c^dlV^— A 



(15) 



We have decided to express the Fokker-Planck equation in terms of the four- velocity 
instead of the four-momentum in order to facilitate the implementation of some 
results used by Rezania and Sobouti in the collisionless case pp. Now f(x^, U l ) is a 
phase density in the six dimensional (x l , U l )-space and a is the Newtonian differential 
cross section. The l.h.s. of (TT2j) in 1PN approximation can be written as p] 
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where U°/c, determined from fllip . is given by 
U° v 2 $ 



v 
2^2 



In consequence, equation ffl2l can be put in the form 
(£jv + £pn)I = A 



.9 1 d 2 

(M l ) + o^H777(/ S 



(16) 



(17) 



(18) 



dU l w ' 2 dU*dW 
where Cn and Cpn are the classical and post-Newtonian Liouville's operators, given 
by 
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(19) 



id® 

3 M ^ 

v + m 



_d_ 
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and A = c/U°, up to order c 2 , equals to 

v 2 $ 



A = 1 



2c 2 



(20) 



(21) 



According to ( fl8l) . we need an expansion of the Fokker-Planck operator (the term in 
parenthesis in the r.h.s.) up to order c~ 2 . Here we have to take into account that 
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Q~ X v 3 {goi + 9ikV k /c) fori^k 



dv 3 
dU 1 

where 



-Q-'lc'U ' 2 +J2 k ^v k (g ok + g lk v l /c)] for i = k 



Q = U°(g 00 + g ol v l /c). (22) 
Then we obtain 

d dv 3 d 



dV i dU* dv 3 

with 



3v 2 $ 

^ 1 -2? + ?' < 24 > 
Now we shall obtain A 1 and S 1 - 7 up to order c~ 2 . It is useful to consider that 
in the center-of-mass system U® = U°, JJ\ = —U l and, as a consequence of the 
condition fTTTj) . we have that g 00 = (700* , 9oi — ~9oi*, 9ij — 9tj*- This will help us to 
simplify the calculation of quantities like AUl/Uo*, d 3 U* and U*U a . In fact, after 
some calculations, we find 

AUi ( 1 2$\ Avl , 1 2$ j3rr / 5v 2 3$\ ,, 



C/*o \ c 2 J c 1 c 2 ' \ 2c 2 c 2 

and 

V^o? - = c (^1 + — - — J |v-v*| 
By introducing the above relations in (f!ll) - (JT5|) . we obtain up to order c~ 2 

A i = J 7 |v-v»|/»Au:a^M\v, (25) 

B 13 = J r]\v-v^AvlAviadQd 3 v„ (26) 



with 



, 7v 2 11$ , 4v 2 12$ 



From equations (j2"5 ]) - (l26"]) one can see the relation between A\ B %3 and the classical 
Rosenbluth coefficients. In order to do this one has to assume that /* can be replaced 
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by a distribution function of a homogeneous field stars population in equilibrium and 
to consider the differential cross section corresponding to the gravitational inverse- 
square law of force. That is 

/(x,v„t) -> *(v,), b = G 2 m 2 |v-v*|- 4 sin- 4 (#/2), (28) 

where ^(v*) = ^* is the field star DF and 9 is the scattering angle as measured in 
the center-of-mass frame. By implementing (1281) in (T2"5]) - (l2"6]) one can see that the 
c°-factors of A 1 and B^ are equivalent to the Rosenbluth's coefficients A l N and B 1 ^, 
given by ([SI)-© (details of this calculation can be seen in j5]). Therefore, we can 
write 

A* = -A* N - A' PN , B» = -B% - B$ N , (29) 
where A l N and B 1 ^ are given by ([2]) and 

i G 2 m 2 f f7v 2 _ y{v*)Avidnd 3 v* 
ApN ~ c 2 J \2c* c 2 J |v-v,|3sin 4 (0/2)' W 

l3 G 2 m 2 f f4v 2 12$ \ V(v,)AviAvidnd 3 v. 

PN c 2 J \ c 2 c 2 J |v-v,|3 s in 4 (0/2) ' 1 j 



Finally, introducing ( 1291 and ( 1231) in ( 1181) . we find that the 1PN approximation of 
(1121) can be written as 

( C N + C PN)! ^-^ N f) + \^(B%l) 

+ j2 [^.(4,/) + K-a(B%f)} , (32) 
where we have defined the operators J7i and /Qj as 

JJ =(, 2 -2$)A + ^(2 + ^, (33) 

. 2 <9 2 3vV <9 2 

Kii = [v 2 — ] — — + 



dvidvj 4 dvidvj 

+ ^( 2+ «W + — s?- (34) 
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The relation (|32|) . in contrast with (jlj), is not a manifestly covariant kinetic equation 
that can be interpreted as an 1PN extension of the classical FPR equation (OQ). 
There, the Newtonian and post-Newtonian contributions appear separated both in 
the l.h.s and in the r.h.s. We note that the post-Newtonian fields ip and & contribute 
only through Cpn and do not appear in the collision term. In the next sections we 
will see that a derivation of the kinetic equation from microscopic dynamics reveals 
additional contributions to the diffusion coefficients. The reason is that the 1PN 
interaction has ternary contributions, coming from ip and that do not disappear 
when the BBGKY sequence is closed. 

3. Derivation from the Microscopic Dynamics: The First 1PN-BBGKY 
Equation 

In this section, we will derive the statistical picture of the self-gravitating system 
starting from its microscopic dynamics. In particular, we shall deal with the BBGKY 
hierarchy that formally equals to the Liouville's equation. In particular, we shall 
consider only the first equation of the sequence, from which we will obtain the kinetic 
equation. Since we start from the non-covariant 1PN law of force, the corresponding 
evolution equation will be also expressed in a non-manifestly covariant fashion. As 
it was pointed in the previous section, this is a feature characterizing the 1PN 
approximation. 

Let a system composed by N identical point particles, with mass m, interacting 
gravitationally and moving with velocity <C c. According to the 1PN approximation, 
each particle experiments an acceleration given by 



where gx and Tr are the newtonian and post Newtonian gravitational forces (per 
unit mass), respectively (the subscript K, motivated by the subsequent use of 
Klimontovich functions, indicates that we deal with exact fields instead of mean 
fields). They are given by [22J 




(35) 



N 




(36) 



^-|-(2<^)-§ + vx(|-x f ) 
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(37) 



+ 3v 



dt 



4v 



d_ 
<9x 



,0$ 
<9x 



The last equation represents the relativistic correction to the force in the 1PN 
approximation, which includes the contribution of the post-Newtonian potentials 
£ and ip . In this case (identical point-like masses), we can obtain an explicit form 
for Fx, in terms of positions and velocities. We have found that this post Newtonian 



contribution can be written in a very suggestive form (see Appendix A): 



N 



K 



5>( 



N 

w, Wi) + 2J 2J "*"( x ' x 

i=l j^i 



(3? 



where w = (x, v) (this notation will be used henceforward). Here, the total 
relativistic force is shown as the result of two contributions: a velocity-dependent 
binary interaction term A, and a velocity-independent ternary interaction term T. 
They are defined as follows 

~4mG 



A(w, w,: 



Gm 



+ - 



+ v 



2v; + 4v ■ Vi 



(4v - 3v<) • r, 



(39) 



TT(x, Xj, Xj 



GV 
2c 2 r T . T 



-7- 



with 



,(40) 



(41) 



As we shall show, A leads to the incorporation of the two-point correlation function 
in the statistical description. On the other hand, T causes the apparition of third 
order correlation patterns, which is an essential difference with the purely Newtonian 
case. 

In order to introduce the statistical description of the evolution for this system, 
we use the KD approach [201 EI]- Such formulation of non-equilibrium statistical 
mechanics, had been widely used in the probabilistic treatment of systems dominated 
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by interactions of kind (1351) (see [23j [2U [25], as examples). The method starts by 
introducing the microscopic one particle phase space density (KD function): 

TV 

Mx, v,f) = £ 5[x - Xi(t)]«y[v - Vi (*)], (42) 

i=l 

where 5 is the 3-dimensional Dirac delta function. Clearly, this function satisfies the 
following evolution equation (KD equation): 

dfx . df K df K df K 

^ r + v-^ r + gx-^ r + r,.^- = 0, (43) 
where gx and T K are given by fl36|) and fl38l) . or by the equivalent relations 

g*(x,t) = -Gm j d 6 w'f K (w',t) r -^, (44) 

rjf(w,*)= J d 6 w'f K (w' } t)\(w,w') 



(45) 



+ y d 6 w'rf 6 w'7^(w',t)/ i ,(w",t)T(x,x',x"). 

The one particle distribution function (representing the system state), as much 
as the first BBGKY equation (describing its temporal evolution), are obtained 
averaging the KD function and the KD equation, respectively. Indeed, it is easy 
to see that the average of fx is 

(/*r(w,f))=/(w,t), (46) 

where / is the usual one particle distribution function, which represents the behavior 
of an average particle in the system. Also, there exist relations connecting correlation 
functions, / and f K |H1 [23] : 

UkS'k) = ff + <?2(w, w', t) + 5(w' - w)f, (47) 
(fKfxfx) = fff" + ^(w, w', w", t) + fg 2 (w', w", t) + f'g 2 (w, w", t) 

+ f"g 2 (w, w', t) + S(w - w') [ff" + <fc(w, w", t)} 

+ 5(w - w") [ff + g 2 (yr, w', t)] + 5(w - w')5(w - w")f 



+ 5(w'-w")[ff + g 2 (w,w',t)}. 



(48) 



Fokker-Planck-Rosenbluth-Type Equations in 1PN Approximation 



12 



Here g 2 and g% are the two and three point correlation functions, respectively, and 
/ = /(w,t),/ / = /(w',t),etc. 

As it is usual, we assume that / must satisfy a kinetic equation describing its 
evolution. In order to derive such equation, we start from the BBGKY first equation, 
which is obtained by taking the average of (1431) . using relations (|46l) - fj48l) : 

9f v 9f 9/ df d 

dt <9x ^ <9v dv dv ' 

where g, T (newtonian and post Newtonian average gravitational fields) are given 
by 

g(x,t) = -Gm f dV/'^fS (50) 

r(w,t) = J d 6 w7 / A(w,w') + y d 6 w'd 6 w"ff"r(xL,x',x"), (51) 
while in the r.h.s. (collision term), F has the form 



A(w,w') — Gm-^- 

^ xx' 



J d e w'g 2 (w,w',t) 
+ J d 6 w'd 6 w" [g 2 (w',w",t)f + g 2 {w,w",t)f 



(52) 



+^ 2 (w, w', t)f" + ^ 3 (w, w', w", t)\ T(x, x', x"). 

Equation (|49|) is the first equation of the 1PN-BBGKY hierarchy. The l.h.s., 
(convective term) depends on the mean fields g, T, and has the same structure 
of a mean field kinetic term (the l.h.s. of a Vlasov equation [26], [27]). More over, 
it is easy to see that this convective term is equivalent to the one characterizing eq. 
(|18|) . However, the collision term is given in a non-standard form. We can see it 
as a sum of two contributions: (i) a two-point term, characterized by g 2 and binary 
interactions g, A; (ii) a three-point term, characterized by the ternary interaction 
T and ternary correlation patterns of the form fg 2 and g%. This last contribution 
makes non trivial the task of close the hierarchy, in order to derive a tractable kinetic 
equation. However, as it will be shown later, we find that is possible to do this, 
establishing some reasonable assumptions about the correlation functions behavior. 
In particular, we will assume that the star cluster relaxation process follows a sort of 
WCG correlation dynamics. Such assumption comes from the fact that the classical 
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FPR equation can be derived directly from the first BBGKY equation by introducing 
a two point correlation function corresponding to a WCG in hydrodynamical regime. 
The validity of such statement will be the subject of the next section. 

3.1. An Alternative Derivation of the Classical FPR equation: the 
Hydrodynamical- WCG Approximation 

As it was pointed out in sectjl], weak encounters play a dominant role in the star 
cluster relaxation process [SI El- This fact suggests that we can expect a close analogy 
between the star cluster relaxation process (dominated by weak encounters) and the 
relaxation process of a WCG (based on weak interactions). In fact, we can show 
that the introduction of a WCG correlation function (in the hydrodynamic regime) 
in the collision term of the classical first BBGKY relation, leads to the usual FPR 
equation. Before to do this, we expose briefly some relevant aspects of the WCG 
correlation dynamics. 

A classical WCG is commonly described as a set of particles interacting via 
gaussian potentials of the form [13] 

V(r) =V exp[-(r/a) 2 ], (53) 

where r is the inter-particle separation, V a represents the potential's maximum value 
and a is an interaction length (for r greater than a the interaction practically 
vanishes). The evolution in time, for this class of systems, is described by the proper 
kinetic equation to order two (see (Hj, p. 574) and there is an aspect we are specially 
interested, around such equation: the explicit form of correlation functions. Here we 
note that (to order two) g%, etc, vanish and that g% is given by the relation [Tj 



f°° OV(w w' t) 
<? 2 (w, w', t) = J dr 1 ^ ' ] • (V^ + rV«0 [//'], (54) 

where we have used the notation 


Vaw* = -k q-;, V„„- e — u w /=v-v. (55) 

OX OX' (TV OV' 

and 

V{w,w',t) = V e- {lr ^~ TU ^ l/a)2 . (56) 

Moreover, if we consider that the WCG is in the hydrodynamic regime, then the 
spatial variation of / is characterized by a hydrodynamical length Lh, defined as 



/ 

L h = max > a. (57) 

\0f/0x\ 
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Since L h ^> a, the system can be considered practically homogeneous (/ is spatially 
independent), over a distance of the order of a. We may take advantage of this fact 
in order to simplify (l54j) . In the Taylor expansion of /' at x' around x, 

/(x', v\ t) = /(x, v', t) - r xx , ■ df{ *£' ,t] + ■ • • , (58) 

the second term is of the order a/Lh compared to the first (the next terms are 
smaller). So we can take /(x', v', t) « /(x, v', t) and then d[ff']/dx' w 0. Performing 
a similar procedure, starting from the expansion of / at x around x', we also 
obtain d[ff']/dx ~ 0. This means that, in the hydrodynamical limit, we can set 
V xx i[ff] ~ and replace /' by f(x,v',t) in the remaining factor V OT ' [//'] of eq. 
(1541) . Thus, by carrying out the integral (154j) . #2 can be written as 

<fe(w, w',t) = 0(w, w') ■ V w /[/(x, v',t)/(x, v,t)], (59) 

with 

£(w, w') = ^e-^l«? W + ( g wUw , 

I V 7r 

-u^r^/a) [1 + Erf(Qw)] e Q »»' } , (60) 
where Erf(Qw) is the error function and 

Qww' = • (61) 

Now let us contemplate the implications concerning with the implementation of 
(|59~I) in the collision term of the classical first BBGKY equation, which we write as 

-^-■F N = -Gm [ d 6 w'£ 2 (w, w', f)£f* (62) 
dv J r% x , 

Again the subscript N denotes that we are dealing with the Newtonian case. In 



Appendix B| we show that after the introduction of (15911 in the above equation, it 
reduces to 

= (63) 

where we have introduced the definitions 

A i N = - 4^ 2 aV GmJL [ d 3 V f ^ ,V,t ) , (64) 

OVi J |v — v'| 

B%= -271 3 / 2 aV Gm^- J dV/(x, v', t)\v - v'|. (65) 
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From the above relations, we immediately note that (|63|) equals to the classical FPR 
collision term (JTJ by choosing V Q as 



V = 7= In A, 

av7r 



and assume that 

/(x,v',0~*(v'), (67) 

where ^/(v'), in the context of the FPR equation derived in §8.1 of [1], is a field 
star distribution function. With fl66|) - fl67|) . A l N and are exactly the Rosenbluth 
diffusion coefficients (EI)-([3j) , describing the drift and diffusivity generated by a 
homogeneous sea of background particles of mass m and distribution function \l/(v'), 
over a test star with the same mass and distribution function f(w,t). 

The assumption (!6"TI) is provided by the hydrodynamical approximation, which 
establishes that the system is practically space independent over the correlation range 
a. Moreover, as it is usual, we demand that the homogeneous distribution \1/ must 
satisfy the relation 

dV*(vO = (68) 

where V is the system's volume. 

All the above considerations lead us to establish that, in order to obtain the 
FPR equation starting from BBGKY hierarchy, we can choose g% = g± = ■ ■ ■ = 
and a two-point correlation function of the form 

g 2 (w,w',t) = G(w,w') ■ V vv ,[*(v')f(w,t)], (69) 

where Q is given by the equation (160]) . defining previously the constant V Q through 
( |66f ). We have to point out that this model of correlation function, in the context 
of a BBGKY scheme, is in close concordance with models based on scattering-cross- 
section considerations as well as approaches stated on a master equation (stochastic 
pictures). 

3.2. The Fokker-Planck-Rusenbluth Equation in 1PN Approximation 

At this stage we can derive an alternative closed kinetic equation for the self- 
gravitating system in 1PN approximation. Here, an important fact, proceeding 
from our experience of sect. [21 will play a crucial role: since in the 1PN approach 



Fokker-Planck-Rosenbluth-Type Equations in 1PN Approximation 



16 



the momentum conservation for a system of point particles is satisfied if and only 
if each particle obeys the Newtonian equation of motion, the scattering process 
is characterized by a differential cross section associated to the Newtonian law of 
force. In other words: the statistical correlation between two colliding particles, 
at 1PN approximation level, must be described by a two point correlation function 
corresponding to the Newtonian scheme. This means that, as a consequence of 
the statements established in the previous section, we can model satisfactorily the 
relaxation process of the self-gravitating system by introducing (169]) in the collision 
term (l52i) . where we also have to set #3 = 0. 

Then, by introducing @2) in Q33J and setting /' « *(v'), /" ~ #(v") 
(assumption ( IBTj) ) in the collision term (1521 . we obtain (see Appendix C ) 

d 2 



df, df 
at ax 



<9v dv 



dvi ^ 2 dvidvj 



(70) 



where A 1 and B 13 are the lPN-drift vector and the lPN-diffusion tensor, respectively. 
They are given by the relations 



.4' 



8nG 2 m 2 In A 



d_ 



, *(V) 



+ / dVf(v' 



_d_ 

dv) 



_d_ 

dv-: 



(71) 



+ J d 3 vW*(v')*(v") 











dv" dv) 



x,v',v") 



where 



B 



>j 



AnG 2 m 2 In A 



2 



dvidvj 



dVtffv')lv- v' 



2 J M(v')n^(w,v') 



O 13 ( 



d 3 x'A i (w,w , )^ i (w,w / ) 
N 



(72) 



(73) 



+ 2— / dVdYT^x^^x'^V.w') 
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$^(x, v', v") = J d 3 xWT l 5 (x, x', x")^'(w', w"). (74) 

The first term on the r.h.s. of flTTj) and (1721) corresponds to the Newtonian 
contribution (classical Rosenbluth coefficients) and the remainder is the post- 
Newtonian contribution that can be calculated explicitly once we choose a particular 
expression for the field star distribution ^ . In this picture, the post-Newtonian 
diffusion coefficients look like much more involved than in the formulation showed in 
sect. El In this subject, we have to point out that they do not necessarily coincide, 
even vanishing ternary contributions. The form in which A l PN and Bp N determines 
the collision term in the first formulation, differs from the present one. 

4. Concluding Remarks 

We have obtained two formulations of FPR-type equations that can be used to 
model the evolution, in a diffusion approximation, of stellar systems where the 
relativistic effects play a significant role. In the first formulation the contribution 
of third order correlation patterns to the diffusion coefficients is vanished, while in 
the second one this fact is took into account. Since the l.h.s of ( 1521 and ( 170|) are 
equivalent, one can expect they be consistent with the post- Newtonian equations of 
hydrodynamics in general relativity pQ. However the generalization of the Eulerian 
equations of Newtonian hydrodynamics, derived from each approach, will differ in 
the term associated to a drift acceleration, since it will be determined by the diffusion 
coefficients. 

Since both (|32|) and (!70|) are expressed in a non-manifest ly covariant fashion (as a 
consequence of the 1PN scheme) its physical interpretation, as well as their differences 
with the usual FPR equation, can be elucidated from a classical Newtonian point 
of view. We have to point out that they have a common point: Once we choose an 
adequate field star distribution function, the diffusion coefficients (f30 j) -( [3Ti) or (JTTJ)- 
(!72j) can be calculated explicitly, making possible to find numerical solutions of (!32|) 
and ( 1701) (for example, using the formalism showed in [28J). 

The KD formalism allows us to derive, in a relatively simple way, the first 
1PN-BBGKY relation, eq. ( )49l) . It constitutes an interesting example about 
the importance of this method, when we deal with the statistical description of 
dynamical systems described by velocity-dependent interactions (see [21]). Perhaps 
the most relevant fact, that can be noticed through (l4T)j) . is the apparition of three 
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order correlation patterns in the collision term. This feature, absent in the purely 
newtonian case, is a consequence of the ternary interaction T (see eq. (1401) ). One 
could go far away and think that more accurate post Newtonian approximations 
(2PN, 3PN and so on) would lead to descriptions based on correlations of greater 
order. It would be not a surprising fact that the Einstein's gravity theory increases 
the correlation level, in a non-equilibrium statistical mechanics description. 

Assumptions fl69|) and fl67j) . reduce f|49l) (in the purely Newtonian case) to 
the usual FPR equation, ([1]). This fact suggests that the local approximation, 
employed in the stochastic approaches of kinetic theory (see for example §8.3 of 
[1]), is equivalent to the hydrodynamic-WCG approximation, used here. Through 
such approximation we model the star cluster relaxation process by the correlation 
dynamics corresponding to a gas dominated by attractive short range interactions, 
V(r) = —(2GmlnA/ay/n)ex i p(—r 2 /a 2 ). The interaction range a, introduced here, 
defines a characteristic scale of distances over which (i) encounters play a dominant 
role, and (ii) the distribution function can be considered homogeneous. One can 
estimate the magnitude of such length, taking into account the considerations used 
in the local approximation, and say that a must be of the order of the minimum 
impact parameter for which such approximation holds. That is, a ~ R/N, where 
R is the system's characteristic radius [4]. According to these statements, the 



equation ( TfOl) . showed in sect. 13.21 describes the evolution of the distribution function 
corresponding to a typical test particle interacting with an homogeneous sea of test 
stars in thermodynamic equilibrium, taking into account 1PN corrections. 

Appendix A. Post Newtonian potentials for N point-like particles 

In the 1PN approximation, $, £ and ip are given by [22] 





o 





V(x,t) 



/ 



r i <9 2 $(x',t) 



+g r 



2 



(x',t) + Gr 



2 



via 



x\t) 



x-x'l 4vr dt 2 



(A.2) 



i 




(A.3) 
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(We denote a, b = 1, 2, 3). For a system composed by N identical point-like particles 
the energy-momentum tensor components, at this order, are 

N 

T 00 (x',t) = m^5(x'-x l (t)), 



i=l 
N 



T 00 (x',t) =mJ2 
i=i L 

l * 
T 0a (x', t) = m - Mt)) 



8(x' - Xi(t)), 



8=1 

N 



T ab (x',t) = m£>?r#(x' -*(*)). 



i=i 



Introducing (IA.4I) in flA.ll) . and flA.61) in flA. 3[) . we obtain 



iV 



$(x,t) = -^ 



Gm 



|x-Xi(t)|' 



N 



AGmvi 
i |x-Xi(t)| 



(A.4) 
(A.5) 
(A.6) 
(A.7) 

(A.8) 
(A.9) 



while, introducing flA.5j) and flA.71) in flA.21) (expression T aa indicates sum over a), ip 
takes the form 



G 2 m 2 



+ 



i=i 

Gm <9 2 



x — x 

N 



i=i jyi 
d 3 x' 



^ -^2 | | -^2 -^J | 



4tt ^ 2 J |x-x'||x'-x i (t)|' 
In order to solve the integral on the right hand we introduce the identity [29 

1 1 f d 3 k 



(A.10) 



ik-(x-x') 



that permits to write 



Ix-x'l 2vr 2 / k 2 



1 f d 3 kd 3 k' 



|x - x'llx' - x,-| An 4 k 2 k' 2 



e i(k-x-k'- Xl ) e i(k'-k)-x' 
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The integral over x' of this last expression, in virtue of the identity 

5(k'-k) = -i- / d 3 xV( k '- k > x ', 
(2tt) j J 

equals to 

f/V 2 f(l3k ^<^). (A.ll) 



7T 



k d 

Defining cartesian coordinates, such that k z is in the direction of 

r = x-x i; (A. 12) 

and then, changing to spherical polar coordinates and introducing the transformation 
u = kr = k\x — Xj|, the right hand side integral of flA. Ill) takes the form 

/^ e *-^) = 47r|x-x i | Hdu^. (A.13) 
J k Jo u 

Since 

/■oc 

du 



x sin u 7r 



u 



3 4 ' 



we find that ( lA.llj) reduces to 

= -27r|x-Xi|. (A.14) 



rf 3 x' 



Then, we have to take into account the term 

— |x - x (t)\ - V * _ l v » - ( x ~ x »)] 2 _ x-Xj _ (hn 
dt 2 ' |x — Xj| |x — Xj| 3 |x — Xj| dt 

in (lA.lOj) . Here, in agreement with the order of approximation, we must to take 

= -GmY * ~ 
dt ^— ' |x,- — x 

Finally, we can write ( lA.lOj) as 



V>(x,t) = -G 2 m 2 ££ 

i=l iVi ^ 



A' 



f 



2 1 5C ^^'i 1 1 ~^^~3 

2v; [ Vl -(x-x,)] 2 



By introducing relations ( 1A.8I) . ( 1A.9I) and ( 1A.15I) in ( 1371) . we can write it as (EE 



(A.15) 
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Appendix B. Derivation of relations ( 163 j) - (1651) 

Introducing (|59|) in (|62|) . it can be cast as 

- — ■F N = I dVV^.n r V ra( [/(x,v',t)/(x,v^)], (B.l) 

where the subscript iV indicates the newtonian case and fijv is a second rank tensor, 
defined as 



N 



dr / d 3 r 



<9$(r) 
<9x 



<9\^(r — ru„ 
<9x 



(B.2) 



where we have changed the integration domain x' by r = x — x', and called 
$(r) = —Gm/r. It is convenient to express $(r) and V(r) in the Fourier expansion: 



with 



d 3 k$ fc e ik ' r , 



V(r) 



:kr 



Gm 



2tt 2 A; 2 ' V * 
Introducing (IB. 31) in (1B.2j) . we find 



a 3 K 



%-(W2) 2 



57T 



3/2 



ft 



N 



in 



dr / d 3 ke irk u -'$i.Vi.kk 



(B.3) 



(B.4) 



(B.5) 



The integral with respect to r can be evaluated using the representation (see 
appendix 2 of |30j) 



dre 



7l5(x) ±iV[- 
X 



(B.6) 



where V(l/x) denotes the principal part. Taking into account that S(x) is an even 
function, V(l/x) is odd, and $fcVfckk is even in the vector k, we obtain 



Q N = 8n 4 J d 3 kS(k ■ u vv ,)$ k V k kk. 



(B.7) 



The calculation of the above integral is simplified using spherical coordinates and 
choosing the Z-axis in the u vv > direction (details of this transformation are shown in 
J, sect 11.6). We find that 



ft 



N 



c 



2 T 

U,,,,/± U,„/U„,/ 



;b.s) 
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C = 8vr 5 / dkk 3 $ k V k = -ir 3/2 GmaV . 
Jo 
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(B.9) 



The relation ( IB. 81) allows us to write the velocity divergence of F^v in the Fokker- 
Planck form: 



d F 

(TV 



where we have used the sum convention (i,j = 1, 2, 3) and 

d d ~ 



A 



A' 



rfV/(x,v') 



ry-3 



B% = 2j dV/(x,v')n«. 
Finally, introducing (1F3. 8|) and flB.91) in the last relations, we obtain fl64|) - fl65|) 



(B.10) 

(B.H) 
(B.12) 



Appendix C. Derivation of the post Newtonian diffusion coefficients 

It is possible to write the remaining part of the collision term in a Fokker-Planck 
form, similar to relation ( |63l) . We first write F = Ftv + Fpat, considering only Fpn 
(the post Newtonian part of F). From fl52|) we note that a simplification is introduced 
if we define symmetric and antisymmetric (with respect to x', x") fields 



Y 5 (x,x',x") = -[T(x,x',x") + T(x,x",x')], 
Y A (x,x',x") = i[T(x,x',x") -T(x,x",x')]. 



(CI) 



(C.2) 



With the help of these definitions, remembering that (72 is symmetric and neglecting 
g 3 , F PN can be cast as 

FpAr = J c? 6 w'A(w, w')(72(w, w') 

+ J d 6 w'd 6 w"Y s (x, x', x") [g 2 (w', w")f + 2^ 2 (w, w')f"} . (C.3) 
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Introducing (1691) in this relation, taking into account the assumption (157)1 . by which 
one can set /' ~ \&(v'), /" w \&(v") in the collision term, we find that the velocity 
divergence of the resulting expression can be written as 



Here we have introduced the post Newtonian diffusion coefficients: 

d d 



where 



A 



PN 



dv'j dvj 



B %3 

PN 



n 13 



+ J d\'d\"^(y')^(y") 
-2 [ dV^v')^ 



9 9 



PN> 



d 3 x'A 4 (w,w')£ J (w,w') 



N 



+ 2y I d 3 x'd 3 x"T 5 (x,x',x")^(w,w'), 



(C.4) 



(C.5) 



(C.6) 
(C.7) 



d 3 xWT s (x, x', x")^'(w', w"). (C.8) 
We obtain (170!) adding (1C.4I) and (IB.lOh (choosing V Q according to (I66j)). 
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